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ABSTRACT 

In a recent paper jBrown & BattvdbOl II) . we proposed the use of integrated polarization measurements 
of background galaxies in radio weak gravitational lensing surveys and investigated the potential impact on 
the statistical measurement of cosmic shear. Here we extend this idea to reconstruct maps of the projected 
dark matter distribution, or lensing convergence field. The addition of polarization can, in principle, greatly 
reduce shape noise due to the intrinsic dispersion in galaxy ellipticities. We show that maps reconstructed 
using this technique in the radio band can be competitive with those derived using standard lensing techniques 
which make use of many more galaxies. In addition, since the reconstruction noise is uncorrelated between 
these standard techniques and the polarization technique, their comparison can serve as a powerful check for 
systematics and their combination can reduce noise further. We examine the convergence reconstruction which 
could be achieved with two forthcoming facilities: (i) a deep survey, covering 1.75 deg^ using the e-MERLIN 
instrument currently being commissioned in the UK and (ii) the high resolution, deep wide field surveys which 
will eventually be conducted with the Square Kilometre Array (SKA). 

Subject headings: cosmology: theory - dark matter, gravitational lensing: weak, radio continuum: galaxies, 
polarization 



1. INTRODUCTION 

Weak gravitational lensing has become one of the most 
powerful techniques for investigating th e distribution of dark 
matter on cosmological scales (e.g. Ma ssev et al.l 120071 : 
iHevman^e t al. 2008; Schrabbacket al. 2010!). The majority 
of weak lensing studies to date have been conducted in the 
optical bands since large numbers of galaxies are required 
to reduce the shape noise associated with the intrinsic ellip- 
ticities of galaxies. Observationally, the pace of progress in 
the optical is set to continue with forthcoming wide-field sur- 
veys using a number of purpose-built instruments, e.g. DE^ 
KID^i Euclicfl WFIRS113 and Pan-STARRSQ. 

As the statistical precision of optical surveys improves, 
greater control over potential systematics is required. Ma- 
jor systematics of concern include instrumental effects such 
as anisotropics in the point-spread function (PSF) of the tele- 
scope and astroph ysical systemati cs suc h as intrinsic galaxy 
ahgnments (e.g ICrittenden et al.ll200lHCatelan et al.ll200Ti: 
IHirata & Seliakll2004D . 

In addition to a new generation of wide-field optical sur- 
veys, a suite of powerful next-generation radio instruments 
are due to come online in the near future. Radio instruments 
offer a number of potential advantages over optical surveys 
for lensing studies. In particular, radio interferometers do not 
suffer from complicated PSF effects, their synthesized beams 
being precisely and solely determined by the telescope array 
configuration and direction of observation. 

A further unique advantage offered by measuring lensing in 
the radio band is the polarization information which is usually 
measured in addition to the total intensity in radio surveys. 



Previous authors have exploited the fact that the polarization 
position angle is unaffected by lensing in order to measure 
gravitational lensing of distant quasars (Kronberg et al. 199U 
1996: .Burns et al..,2004.) . In a recent paper ( Brown & Battyg 
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201 II) . we showed how one could extend this idea to measure 
cosmic shear. The technique relies on there existing a rea- 
sonably tight relationship between the orientation of the inte- 
grated polarized emission and the intrinsic morphological ori- 
entation of the galaxy. The existence of this relationship needs 
to be established for the high-redshift star-forming galaxies 
which are expected to dominate the radio sky at the /iJy flux 
sensitivities achievable with forthcoming instruments. How- 
ever, such a relationship certainly exists in the local Universe 
(Stil et al. 2009) and it is reasonable to assume that it persists 
to higher redshift. 

A key difference between the polarization technique and 
standard techniques for measuring lensing is that the former 
does not assume that the ensemble average of the intrinsic 
shapes of galaxies vanishes. It is thus, in principle, able to 
cleanly discriminate between a lensing signal and a possi- 
ble contaminating sign al due to intrinsic galaxy alignments 
(iBrown & Battyeil201 ll) . 

In this letter, we investigate the direct reconstruction of 
the projected dark matter distribution (or lensing convergence 
field) using the radio polarization technique. Our simulations 
are based on the ray-tracing simulations of White (2005). 
We begin in Section |2] by summarizing the salient features 
of these simulations and we describe the parameters we have 
adopted to mimic the e-MERLI>0 (currently being commis- 
sioned with full operation planned for the second half of 20 11 ) 
and SK/^ instruments. In Section [3] we describe how to in- 
clude polarization information in a direct-inversion mass re- 
construction algorithm and we present our simulation results. 
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We finish in Section|4]with a discussion. 



2. SIMULATIONS 



To demonstrate the mass reconstruction technique using po- 
larization information, we have made use of the simulated 
lensing convergence and shear maps of White (.20051^ These 
simulations comprise ss 1000 square degrees of simulated 
sky from high resolution numerical simulations based on a 
ACDM cosmology with parameters, Q,,, = 0.28, Clhh^ = 0.024, 
h = 0.7, (Ts = 0.9 and n., = 1. The simulations consist of 112 
3x3 deg^ fields which are approximately independent. 

The assumed redshift distribution is of the form, n{z) oc 
z^exp(-z/zo)'^ with zo = 1-0 and /3 = 1.5. The median red- 
shift of this distribution is Zm = 1.41. Making use of the 
semi-empir ical SKA Design Studies (SKADS) simulation of 
iWilman et a l. (20080, we have checked that this n(z) agrees 
reasonably well with the redshift distribution of star-forming 
galaxies expected in future deep radio surveys. In particular, 
we measure Zm = 1-02 in the SKADS simulation for a IOct 
detection threshold of 40 /iJy, a sensitivity which could be 
reached over 1 .75 deg^ using e-MERLIN. Large-scale s urveys 
with the SKA might go ~ 10 times deeper (e.g. Bla ke et alj 
l2007h . At this depth, we measure a median redshift of 
Zm = 1-48 from the SKADS simulation. Although the pre- 
dicted median redshift of the proposed e-MERLIN survey is 
somewhat smaller than that of White's simulations, this dis- 
crepancy is sub-dominant when compared to the current very 
large uncertainties in predicting n{z) for /xJy-sensitivity radio 
surveys. 

A further unknown for forthcoming radio surveys is the 
number density of galaxies that will be achieved . Source 
counts at 1.4 GHz are presented in Biggs & I vison' ( 2006) for 
deep fields over relatively small areas of sky. Extrapolating 
from this study suggests a source density of ^ 1 arcmin"^ for 
a detection threshold of 40 /iJy. This source density is also 
consistent with that found in deep VLA + MERLIN obser- 
vations of the Hubble Deep Field (Muxlow et al. 2005). The 
source counts from a more recent study by Owen & Morriso^ 
(j2008i) are significantly higher and suggest a source density 
of 2 arcmin"^ could be achieved. For the e-MERLIN simu- 
lations, we take the average of these two estimates and adopt 
n = 1 .5 arcmin"^ for the number density of galaxies detected 
in total intensity. For the SKA simulations, we assume n = 15 
arcmin"^. 

The polarization properties of the high-redshift star- 
forming galaxies which will dominate these surveys are cur- 
rently unknown. We can only extrapolate from the local Uni- 
verse wh ere the fractional polarization is typically between 
and 20% jStil et al. 2009). These observations of local galax- 
ies also suggest that the polarized emission and the morpho- 
logical orientation are well aligned to an accuracy of between 
and ^ 15 degs. 

For our simulations, we assume that the galaxies are, on av- 
erage, 10% polarized. Noting that a measurement of the ori- 
entation of the polarized emission requires only a third of the 
sensitivity which would be required to measure the shape of 
the emission, we assume further that we can measure a use- 
ful polarization orientation for one third of the galaxies de- 
tected in total intensity. Thus, we assume Mpd = 0.5 arcmin"^ 
for e-MERLIN and «poi = 5.0 arcmin'^ for the SKA. We fur- 

^ Ava ilable at|http : / /mwhite .berkeley ■ edu/Lensing/ThousandT] 
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ther assume that such polarization orientation detections are 
an unbiased tracer of the intrinsic morphological orientation 
of the galaxy with a scatter of arms = 7 degs. (See Section 2 
of Brown & Battve 201 1 for further discussion on the current 
observational status for polarization and possible values for 
these parameters.) 

To mimic the proposed e-MERLIN survey, we extract the 
central 1.75 deg^ from each 3x3 deg^ simulation. An ex- 
ample of the convergence distribution for one field is shown 
in the top left panel of Fig. [T] where a number of large-scale 
structures and filaments are clearly visible. Using the param- 
eters detailed above, we generate background galaxies ran- 
domly distributed across the field of view. Each galaxy is 
assigned an intrinsic ellipticity, e = ei+ie2 where ei and £2 
are drawn from zero-mean Gaussian distributions with width 
fJe = 0.30. For each galaxy, we then add the simulated shear 
to the intrinsic ellipticity, interpolating the shear from the nu- 
merical simulations. 

For the polarization technique, we additionally simulate the 
observed orientation of the polarized emission of each galaxy 
according to a°'°'' = a'"^ + a™'^ where a'"' = (1 /2)tan-'(e2/ei) 
is the intrinsic position angle of the galaxy and a™'^ is the 
random noise in the polarized orientation — intrinsic position 
angle relationship. 

3. MASS RECONSTRUCTION USING POLARIZATION 

To reconstruct the convergence field from the simulated 
galaxy catalogues, we use a variant of the Kais er-Squires 
direct-inversion technique jKaiser & Squireslll993h . Our im- 
plementation is similar to the real-space algorithm described 
in Seitz & Schneider (1995). To perform the reconstruction 
using the total intens ity only (i.e. using the sta ndard lensing 
technique), following lSeitz & Schneiden(ll995 '). we first esti- 
mate the mean shear field on a regular grid covering the entire 
field. We estimate the shear in each grid cell as (the indexes, 
{i,j } denote the grid cells) 



(1) 



where the sum is over all galaxies and we have applied a 
weight of 



Wk = exp 



(A6l)2 



(2) 



to each of the observed ellipticities, e^*'^ We use a smooth- 
ing scale of A9 = 2.0 arcmin. The convergence field is then 
estimated as a convolution. 



k(e) = ^Re[I?(0-0,-,,)7*(0,-.,)] 



where the kernel, V is given by 



] + 2i9 
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Including the polarization information in this estimator is 
straightforward and simply involv es using a weighted vers ion 
of the shear estimator described in lBrown & Battv^ (1201 Ih to 
estimate the shear at each grid point. Explicitly, the shear in 
each grid cell is now estimated as 



7 = A-'b, 



(5) 
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Figure 1. An example reconstruction of the dark matter distribution in a 1.75 deg^ region of the simulations. The input convergence field is shown unaltered in 
the upper left panel and smoothed with a 2 arcmin Gaussian kernel in the lower left panel. The simulated reconstruction using total intensity only for observational 
specifications similar to the proposed e-MERLIN survey (with n= 1.5 arcmin"^ and eims = 0.30) is shown in the upper central panel. Assuming one third of these 
galaxies are sufficiently polarized to yield an intrinsic position angle tracer that is unbiased with a scatter of arms = 7 degs, the reconstruction obtained using 
polarization information is shown in the lower central panel. The reconstructions obtainable with the SKA are shown in the right-hand panels (upper panel — 
using total intensity only; lower panel — with polarization). For the SKA reconstructions, we have assumed n = 15 arcmin"- and that one third of the galaxies 
are useable in polarization with Qnns = 7.0 degs. Integrated convergence measurements for the two structures labelled "Cluster 1" and "Cluster 2" are given in 
Tabled 



where the matrix, A, and the vector, b, are given by 



(6) 



(7) 



('sin2af% 



-cos2af') 



is the direction at 45° to 



Here, 11^ 

the estimate of the intrinsic position angle of the galaxy as 
provided by the polarization and is a weight. 

Using a weighting scheme in equation Q is now crucial - 
the small galaxy number density we are dealing with means 
that "shot noise" is a severe problem for the polarization es- 
timatorF^ Had we used an unweighted shear estimator (i.e. 
simple binning) then a significant fraction of our grid cells 
would contain no galaxies for the polarization estimator and 
this would present problems for the convolution operation. By 
using weights which are non-negligible at the typical separa- 
tion of polarized galaxies, we can recover the shear at each 
point in the grid, albeit at the expense of a smoothing of the 

We discriminate here between "shape noise" which is the dispersion in 
the intrinsic shapes of galaxies and "shot noise" which is the noise introduced 
due to the finite number of galaxies and hence the sparse sampling of the 
underlying shear field. 



shear field. In practice, we use the same weighting scheme 
as for the standard case (equation |2]l with the same smooth- 
ing scale. Ad = 2.0 arcmin, in all cases. Once the shear field 
has been estimated, the convergence, k(6) is found as in the 
standard case (equation[3]). 

An example of the mass reconstruction performance for 
both the standard estimator and the polarization estimator is 
shown in Fig. [T] The simulated e-MERLIN reconstructions 
have been performed on a 64 x 64 grid of 1.25 arcmin pix- 
els. For the much higher fidelity SKA simulations, we use a 
128 X 128 grid of 0.625 arcmin pixels. We see that the ma- 
jor mass concentrations can be successfully identified in the 
simulated e-MERLIN reconstructions while the SKA maps 
also trace the filamentary structure of the numerical simula- 
tions extremely well. For the observational parameters we 
have adopted, the reconstructions obtained using the polar- 
ization estimator are of similar sensitivity to those obtained 
using the standard estimator. In Table [T] we list the recon- 
structed convergence field integrated in apertures of radius 
3 arcmin for the two largest structures in the field, centred 
on {x,y) = (32.6,42.8) arcmin and {x,y) = (56.5,24.3) arcmin. 
The errors are determined by repeating the analysis on empty 
fields and measuring the root-mean-square convergence field 
integrated in identical apertures. 
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Table 1 

Integrated convergence measurements for the structures 
labelled "Cluster 1" and "Cluster 2" in Fig.[T] 

Simulation Standard estimator Polarization estimator 



e-MERLIN: 
Cluster 1 
Cluster 2 

SKA: 
Cluster 1 
Cluster 2 



2.36±0.19 
1.26±0.19 



2.29 ±0.08 
0.92 ±0.08 



2.31±0.18 
1.30±0.18 



2. 13 ±0.07 
0.98 ±0.07 



We have also measured the shear correlation functions from 
the simulated maps. The correlation functions are defined as 



Ci(0)=(7[(0)7i(0+A0)) 

C2(0)=(72(0)72(^+A0)) 



(8) 



where the angled brackets denote an average over all map- 
pixel pairs separated by and 7J and 72 are the estimated 
shear components rotated to a coordinate system aligned with 
the vector joining the two map-pixels. The cross-correlation, 
CiiO) is expected to vanish in the absence of parity-violating 
effects. 

We have measured the total correlation function, C{6) = 
Ci{6) + C2{0) using both the standard and polarization esti- 
mator and have also combined the results from the two tech- 
niques. Note that the combination of the two estimators is 
powerful since the dominant source of measurement noise 
is different in the two cases — for the standard estimator, 
the shape noise dominates while for the polarization estima- 
tor the noise associated with a^ms dominates meaning that 
there is useful extra information in the cross-correlation of the 
shear fields estimated using the two techniques. The cross- 
correlation can also be useful in mitigating systematics which 
would likely be different for the two techniques. 

The results, averaged over simulations, are plotted in Fig.|2l 
Naively, it looks as though a detection of cosmic shear is pos- 
sible with the proposed 1 .75 deg^ e-MERLIN survey. How- 
ever, the measurements on different scales are strongly corre- 
lated and the formal detection significance is only 1 Aa. Note 
that the dominant contribution to the errors is sample vari- 
ance. Instead of doing cosmology, if we ask the question how 
well can we measure the lensing signal in a typical 1 .75 deg^ 
field, then the significance is much larger 6.6<j). 

4. DISCUSSION 

In t his letter, we have extended the techniques of 
iBrown & Battyd (1201 Ih to the lensing reconstruction of dark 
matter maps from polarized radio surveys. For the parame- 
ters which we have adopted, the technique looks extremely 
promising and could make future radio surveys as powerful 
as optical lensing surveys containing many more galaxies. 

There is, of course, a great deal of uncertainty with respect 
to the values one should adopt for the polarization properties 
of star-forming galaxies at /iJy flux densities. In particular, 
the typical fractional polarization and the degree to which the 
orientation of the polarized emission is aligned with the in- 
trinsic morphological orientation are currently unknown. We 
note also however that there is a trade-off in these two parame- 
ters. For example, one can imagine including a larger number 
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Figure 2. The total shear correlation function measured from the e-MERLIN 
simulations. Thi'ee sets of points are plotted. From left to right they have been 
obtained using the standard shear estimator, using the polarization estimator 
and their combination. The points plotted as open circles show the cross- 
con'elation, C3 (0), which is consistent with zero. Comparison of the en'ors on 
this latter measurement with those on the measurement of the signal gives an 
indication of the relative sizes of the measurement noise and sample variance 
contributions. The errors plotted are those appropriate for a single realization. 
Measurements on different scales are heavily correlated and the significance 
of the detection of cosmic shear is in fact only \Aa. 

of galaxies whose polarized emission and intrinsic orientation 
are poorly aligned (large «poi and large anii,s) or one could be 
very selective about which galaxies one includes, choosing 
high fractional polarization objects whose polarization orien- 
tation closely traces the major axis of the galaxy (low Wpoi and 

low arms). 

In terms o f forecasting constraints for future surveys, 
iKaised (1199 8*) has shown that the errors obtainable on the 
shear power spectrum with a survey covering a fraction of 
sky /sky and with a galaxy number density n are given by 



AQ = 



C, ^rms 

n 



y {21+ l)/sky 

The equivalent expression for the polarization case is 



AQ = 



(21+ l)/s 



sky 



16a 



2 e2 

rms rms 



Wpol 



(9) 
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In Fig. [3] we have plotted these predictions for some repre- 
sentative numbers chosen to mimic large-scale surveys to be 
conducted with the SKA. If the polarization properties of star- 
forming galaxies assumed here turn out to be reasonable, then 
it will be very interesting to apply the polarization techniques 
described here and elsewhere to forthcoming high-resolution 
/ijy-sensitivity radio surveys. 
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Figure 3. Predicted per-multipole uncertainties on the weak lensing power 
spectrum for a 5000 deg^ survey with n= 15 arcmin"^ and emis = 0.3. The 
forecasts using the polarization estimator assume that one third of these 
galaxies are detected sufficiently well in polarization such that arms = 7 degs. 
The cosmic variance limit is also shown for comparison. 



